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ABSTRACT 

It is generally believed that turbulence has a significant impact on the dynamics and evolution 
of molecular clouds and the star formation which occurs within them. Non-ideal magnetohydro- 
dynamic effects arc known to influence the nature of this turbulence. We present the results of a 
suite of 512 3 resolution simulations of the decay of initially super- Alfvenic and supersonic fully 
multifluid MHD turbulence. 

We find that ambipolar diffusion increases the rate of decay of the turbulence while the Hall 
effect has virtually no impact. The decay of the kinetic energy can be fitted as a power-law in 
time and the exponent is found to be —1.34 for fully multifluid MHD turbulence. The power 
spectra of density, velocity and magnetic field are all steepened significantly by the inclusion of 
non-ideal terms. The dominant reason for this steepening is ambipolar diffusion with the Hall 
effect again playing a minimal role except at short length scales where it creates extra structure 
in the magnetic field. 

Interestingly we find that, at least at these resolutions, the majority of the physics of multifluid 
turbulence can be captured by simply introducing fixed (in time and space) resistive terms into 
the induction equation without the need for a full multifluid MHD treatment. 

The velocity dispersion is also examined and, in common with previously published results, 
it is found not to be power-law in nature. 

Subject headings: MHD ISM: kinematics and dynamics ISM: magnetic fields methods: numerical 
turbulence 



1. Introduction 

Turbulence is recognized as a possible source of 
support against gravitational collapse for molec- 
ular clouds. The precise role and source of the 
observed motions interpreted as evidence of tur- 
bulence in these clouds has been studied ex- 
te nsively by many researchers (see the reviews 
of Mac Low fe Klessen 2004 ; Elmegreen fe Scalol 
20041 ). Clearly, if turbulence can support molecu- 
lar clouds then it can influence star formation in 



terms of rate, efficiency and initial mass function 
(|Elmegreenlll993l : iKlein et al.ll2003l ). 

Many studies of turbulence in molecular clouds 
have focused on ideal magnetohydrodynamics 
(MHD) as an app roximation of the phys i cs govern- 



ing t h is system (IMac Low et al.l 119 98; Mac Low 



19991: lOstriker et al. 


2001; Vestuto et al.l 


2003 


Gustaffson et al. 2006 


: Glover & Mac Low 


2007 



Lemaster fc Stone! 1200a l2009h . The assumption 
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of ideal MHD, while desirable for technical rea- 
sons, is perhaps risky in the context of turbulence. 
The reason for this is that while ideal MHD is 
valid in molecular clouds on fairly large length 
scales, on shorter length scales no n-ideal effects 
are thought to become significant (|War diel l2004t 
Oishi fc Mac Lowll2006h . Given that turbulence m 



3 dimensions involves the transfer of energy from 
large scales to ever smaller scales, the assumption 
of ideal MHD will be invalid below some critical 



1 



spatial scale and the correct nature of the energy 
cascade may not be observed at this range. 

The most important of the non-ideal effects for 
molecular cloud dynamics is ambipolar diffusion 
Some authors (lOishi fe Mac Low 2006 ; Li et al 



20081: iKudoh fc Basul 120081) have studied driven 



MHD turbulence in the presence of ambipolar dif- 
fusion. All these authors find that ambipolar dif- 
fusion produces significant differences in the prop- 
erties of the turbulence. 

While most likely of lesser significance, it has 
been suggested that although the Hall resistivity 
is generally at least an order of magnitude lower 
than the ambi polar resistivity in molecular clouds 
(IWardldl200l . its effect should not be ignored. 
Although relatively weak, it is capable of induc- 
ing topological changes in the magnetic field which 
are quite distinct to any influence caused by am- 
bipolar diffusion. In support of this assertion, we 
note that researchers working on reconncction and 
the solar wind have studied the Hall effect in the 
context of turbulence and found that, although 
the overall decay rate appears not to be affected, 
the usual coincidence of the magnetic and veloc- 
ity fiel ds seen in MHD does not occur at small 
scales dMatthaeus et al"1l2003t iMininni et"al1l2006l: 



Servidio et al.l 120071 ). Almost no work has been 
done on comparing the influences of this effect 
coupled with that of ambipo lar diffusion on turbu- 
lence with the exception of iDownes fe O' Sullivan 
([20091 hereafter Paper I). 



In Paper I a series of simulations of decay- 
ing supersonic non-ideal MHD turbulence incor- 
porating both ambipolar diffusion and the Hall ef- 
fect were performed. These simulations, however, 
were constrained in that the resistivities associ- 
ated with each of ambipolar diffusion, the Hall ef- 
fect and the Pederson resistivity were kept fixed in 
both space and time. The authors found that, at 
length scales of 0.2 pc, ambipolar diffusion has a 
significant impact on the decay of the turbulence. 
The Hall effect was less significant in this respect 
but does have an influence on the magnetic field 
at short length scales. Here we present simula- 
tions in which the resistivities are self-consistently 
calculated from the evolution of both the mag- 
netic field and the densities of all of the compo- 
nent species of the fluid. Using these dynamically 
evolving resistivities we study the decay of fully 
multifluid MHD turbulence. This is the first such 



study presented in the literature, with the excep- 
tion of the low resolut i on sim ulations presented by 
IDownes fc O'Sullivanl (|2008h . 

The aim of this work is to examine in detail the 
differences between the decay of ideal MHD turbu- 
lence and that of multifluid MHD turbulence with 
a full tensor resistivity incorporating the effects of 
ambipolar diffusion, the Hall effect and Ohmic re- 
sistivity. We will use the results of Paper I in our 
discussion of these differences as it represents an 
intermediate stage between the calculations pre- 
sented here and those of ideal MHD. This work 
is new in two respects: notwithstanding Paper I, 
no previous work has focused on decaying (i.e. un- 
driven) multifluid MHD turbulence and, in addi- 
tion, no previous work has addressed the issue of 
turbulence in the presence of both ambipolar dif- 
fusion and the Hall effect simultaneously. 

In section[2]we outline the numerical techniques 
used in this work, as well as the initial conditions 
and general set-up for the simulations while in sec- 
tion |3] wc describe the methods used to analyze 
the simulation data. In section |4] we present and 
discuss the results of our simulations of turbulent 
decay. Finally, section [5] contains a summary of 
our results. 

2. Numerical method 



As in Paper I, we use the c ode HYDRA 
(jO'Sullivan fc Downed I2OO6I \2000) to integrate 



the equations of weakly ionized multifluid MHD 
(see section 12 . 1[) . Wc assume that the molecular 
cloud material we are simulating can be treated as 
isothermal and that initially the density and mag- 
netic field are uniform. We use the capabilities of 
HYDRA to extend the physics incorporated in the 
simulations here beyond those presented in Paper 
I so that the turbulence here is fully multifluid 
MHD. 

2.1. Equations and algorithm 

We briefly outline the equations and assump- 
tions in our model h e re but refer the reader to 
O'Sullivan fc Downed (|2006l 120071 ) for a compre- 



hensive description of the underlying assumptions 
for the weakly ionized model of multifluid MHD. 

We assume that the cloud material can be 
treated as weakly ionized. This is clearly valid 
for molecular clouds and allows us to ignore the 
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inert i a of the ch arged species (jCiolek fc Robergd 
20021 : lFallell2003l ). For a system composed of N 
fluids, one of which is neutral, the equations to be 
solved are then 



dfh 

at 



V-( P i 9i )=0, {l<i<N), 
V ■ + a 2 pl) = J x B, 

-V x E', 



dt 
dB 

— + V-( Ql B-B qi ) 



(1) 
(2) 
(3) 



u i p l (E + q- x B) + p l p 1 K il (q 1 - q { ) = (4) 



V B = 0, 
VxB = J, 



(5) 
(6) 



A' 



^2a>ipi = 0, (7) 



where pi, q 1 , a, B and J are the neutral mass den- 
sity, neutral velocity, sound speed, magnetic field 
and current density respectively and 2 < i < N 
unless otherwise noted. .fQi, on and pi (i > 1) 
are the collision coefficients between species i and 
the neutrals, the charged fluid charge-to-mass ra- 
tios and mass densities respectively. Equations ([T]) 
to ([7]) express conservation of mass for each fluid, 
conservation of neutral momentum, the induction 
equation, force-balance for the charged species, 
the inadmissibility of magnetic monopoles, Fara- 
day's law and charge neutrality respectively. 

The electric field in the frame of the fluid, 
E', is calculated from the general ized Ohm's 
law for weakly ionized flu ids (e.g. Falld 2003 ; 
O' Sullivan &: Downed 2006 ) and is given by 



where 



E — Eq + .Eh + Ep 



Eq = (J ■ a )a , 

E H = J x oh, 

E A = -(J x o A ) x a A , 



(8) 



(9) 
(10) 
(11) 



using the definitions ao — /o-B, ah = ./h-B, a a = 
f A B, where / = ^n5/B } / H = rn/B, f A = 
y/rJi/B. ro, ?"H and r A are the Ohmic, Hall and 
ambipolar resistivities respectively and are given 

by 



ro 



1 



(12) 



?'H = 

r A = 



CH 



T 2 



1 ' 



CA 



with the conductivities given by 



TO 



CH 



1 N 



j=2 
JV 



Ok = 



1 x ^ 

~B 1 

i=2 



(13) 
(14) 

(15) 
(16) 

(17) 



where Pi is the Hall parameter for species i and is 
given by 

A = (18) 

As noted bv lFalld (|2003l) and lO'Sullivan fc Downed 
(|2006l ). the main difficulty with standard numer- 
ical techniques for integrating equation ([3]) lies 
with the Hall term. As this term becomes domi- 
nant the stable time-ste p goe s to ze ro. However, 
O'Sullivan & Downes (|2006l l2007h presented a 
novel, explicit numerical method for integrating 
this term such that the limit on the stable time- 
step is not overly restrictive. We use this "Hall 
Diffusion Scheme" in this work. Of course, all 
explicitly differenced diffusion terms give rise to 
a stable time-step which is proportional to Ace 2 , 
where Ax is the resolution of the simulation. To 
ameliorate this we use standard subcycling of the 
Hall terms and super time-steppi ng to accelerate 



the ambipolar diffusion terms (seelAlexiades et al 
19961: lO'Sullivan fc Downesll2006l [2007l) ~ 



Equations (P) - © are solved using a standard 
shock-capturing, second order, finite volume, con- 
servative scheme. The numerical techniques em- 
ployed in this work are slightly different to those 
used in Paper I in one respect: we have altered 
the calculation of the advective fluxe s in equatio n 
© to use the method suggested bv lFalld (f2003h . 
In Paper I these fluxes were derived from interface 
values of the neutral gas velocity and the magnetic 
field. We find that at high resolutions with vari- 
able resistivities the latter method is prone to in- 
troducing grid scale features in the solution while 
the former is not. This undesirable effect was not 
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an issue for the investigations carried out in Pa- 
per I since resistivities were fixed. The downside 
of the described variation between the numerical 
approaches is that it must be considered as a pos- 
sible source of discrepancy in comparisons between 
the results of Paper I and this work. However, in 
order to provide evidence of the small influence, 
we have also run a fixed resistivity simulation for 
this work (see section l2~2j) . 

Equation O is e nforced using the method of 



Dedner et al.l (|2002l) . The effects of the diffusive 
terms in equation ([3]) are then incorporated in an 
operator split fashion. 

2.2. Initial conditions 

We examine the decay of MHD turbulence in 
conditions suitable for dense regions of molecular 
clouds. The conditions we use are similar to those 
used in Paper I. We briefly review them here for 
completeness. 

The computational domain is set up as a cube 
of side L = 0.2 pc. Periodic boundary conditions 
are enforced on all faces of the simulation domain. 
The sound speed is set to 0.55 km s^ 1 , the initial 
density is chosen to be uniform with a value of 
f 6 cm -3 and the magnetic field is also initially 
uniform in the (1, 1, f ) direction with a magnitude 
of 1 mG. For these conditions, suitable conductiv- 
ities are ero = f x 10 10 s~\ oh = IO^s" 1 and 
cr A = lO^s" 1 dWardle k Nell 19991) . We choose 
a 3-fluid set-up for our multifluid simulation: 1 
neutral species and 2 charged species. The den- 
sities, charge-to-mass ratios and collisional coeffi- 
cients of the charged species are chosen in order 
to achieve these conductivities. We choose these 
particular physical conditions with a view to max- 
imizing the influence of the Hall effect in our simu- 
lations (jWardle k Nglll999h . In this way we hope 
to determine whether the Hall effect is ever likely 
to be important in molecular cloud turbulence. 

The initial velocity field is defined to be the sum 
of waves with 64 wave-vectors each with random 
amplitude and phase - i.e. 



64 

E 

3=0 



A a j cos(fcj • x 



(19) 



where a defines the component (x, y or z) of the 
appropriate quantity, A a j and 4> a j are the ran- 
dom amplitudes and phases and x is the posi- 



tion vector. We restrict the velocity field to be 
solcnoidal (i.e. non-comprcssional) . By construc- 
tion the mean velocity over the domain is zero. 

Table [T] presents a complete list of the various 
simulations carried out in this work. The nomen- 
clature we employ in referencing the simulations is 
xx-c where xx denotes the type of physics (e.g. a 
standard molecular cloud run is "mc" , ideal MHD 
is "mhd" etc) and c is the resolution used. The 
initial root-mean-square (rms) of the field is cho- 
sen to be 5 with a corresponding Alfvenic Mach 
number of approximately 1.9. In addition to the 
4 multifluid MHD simulations run at different res- 
olutions, we also run 4 further simulations. The 
first is an ideal MHD simulation (mhd-512) which 
we use for comparison purposes and the other two 
(ambi-512 and hall-512) only incorporate one of 
ambipolar diffusion or the Hall effect, respectively. 
The final case (fr-512) is a fixed resistivity simula- 
tion used to make contact with the simulations of 
Paper I. We use these latter 4 simulations to in- 
vestigate separately the influence of each non-ideal 
effect. 

3. Analysis 

In this section we discuss the method of analysis 
of the output of the simulations described in sec- 
tion l2.2l The aim of this paper is to investigate the 
decay rate of supersonic turbulence in molecular 
clouds. Hence, the main analysis carried out is of 
the kinetic, magnetic and total energy as functions 
of time. These quantities are defined respectively 
as 

e k = / p\q\ 2 dV (20a) 

J domain 

e b = f l*!V-<4>V(20b) 



J domain 

etot = eb + ek 



(20c) 



where V is the volume of the computational do- 
main. 

We also calculate the mass-weighted average 
Mach number, defined by 



1 



M=-{al + a 2 y + a 2 z } 



1/2 



(21) 



where a is the sound speed and the velocity dis- 
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Table 1: Definition of the initial conditions used in the simulations in this work. 



Simulation 


Mach number 3. 


Resolution 


Comment 


mc-64 


5 


64 3 




mc-128 


■5 


128 3 




mc-256 


■5 


256 3 




mc-512 


5 


512 3 




ambi-512 


5 


512 3 


r H =0 


hall-512 


5 


512 3 


r A =0 


mhd-512 


5 


512 3 


Ideal MHD 


fr-512 


5 


512 3 


As mc-512 but for fixed resistivities 



"Initial rms Mach number of the flow 



persions, a a , are defined by 



4. Results 



(Pgg) 
(P) 



1/2 



(22) 



where a is either x, y or z and the angle brackets 
den ote averaging over the co mputational domain 
(see Lemaster h Stonel 2009() . 



In section T4. 31 we present the power spectra for 
the velocity, density and magnetic field for each of 
the 512 3 simulations. These spectra are calculated 
by taking the power spectra in the x, y and z 
directions separately and summing the power over 
the interval k < \k\ < k + Ak (where we take 
Afc = 1). This gives us some insight into the scale 
of structures being formed by the turbulence for 
the various initial conditions and range of physics 
examined. 



Finally, in section 14.41 we calculate the velocity 
dispersion as a function of length scale, I. For 
these purposes we define the velocity dispersion 
to be 

°"(0 = {< ^xCO > domain + < °"y(0 >domain 

domain V (23) 



where 



MO = {(<£>« -<**>?}* (24) 

where < • >i indicates an average taken over a 
cube of side I in the simulation domain and < 
• >domain indicates averaging of the quantity over 
all such non-overlapping cubes within the domain. 



Each of the simulations detailed in table Q] was 
run for one sound crossing time, t c = 3.56 x 10 4 yrs, 
of the simulation domain. All analysis was carried 
out for t > 0.2 t c (i.e. one flow crossing time) at 
which point we expect significant turbulent mixing 
to have taken place and the system's memory of 
the initial state to be largely forgotten. 

As an illustration of the differences between an 
ideal MHD turbulence simulation and a multifluid 
MHD simulation, figure [T] contains plots of the 
density distribution at t — t c in a slice through 
the computational domain for simulations mhd- 
512 and mc-512. It is clear that there is much 
less fine structure in mc-512. Also shown is the 
same slice for simulation fr-512 (i.e. fixed resistiv- 
ities). While the similarities in terms of the levels 
of structure are relatively small between mc-512 
and fr-512, there are clear differences between the 
distributions indicating that calculating the resis- 
tivities self consistently has some impact on the 
dynamics of the system. 

In figure [2] we present plots of the ambipolar 
and Hall resistivities for simulation mc-512 for the 
same times and slices as figure [U It is clear that 
the resistivities vary considerably throughout the 
computational domain with the features strongly 
correlated with the features in the density distri- 
bution. We also show n = -r^i to give an indi- 
cation of the relative importance of each of the 
resistivities. This, as we shall see, is an important 
parameter. Finally, figure [3] is the same as [2] ex- 
cept that the data is taken at time t = 0.2t c - i.e. 
after one flow crossing time. Here we can see that 
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the variation of the resistivities in space is dra- 
matic with, for example, the ambipolar resistivity 
varying by almost 4 orders of magnitude with rj 
varying by around 2 orders of magnitude. 
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Fig. 1. — Log plot of the neutral mass density at 
t = t c in a slice at z = 0.5 for simulations mhd-512 
(upper panel), mc-512 (middle panel) and fr-512 
(lower panel). 



4.1. Resolution study 

Four simulations identical in every way except 
for the resolution were run. Specifically, the reso- 
lutions used were 64 3 , 128 3 , 256 3 , and 512 3 . We 
now focus our attention on how the energy decay 
behaves with resolution. Figured] contains plots of 
the kinetic energy as a function of time for each of 
the simulations in the resolution study. It is clear 
that the lower the resolution, the faster the decay 
- this is what one would expect since lower reso- 
lution results in a higher numerical viscosity and 
hence one expects faster dissipation of energy. 

Simulations mc-256 and mc-512 are, however, 
quite similar in terms of the energy decay with 
a maximum relative difference of around 10% be- 
tween the kinetic energies in the simulations at any 
one time - an almost identical result to that ob- 
tained from the resolution study in Paper I. This is 
notable since in Paper I the resistivities were kept 
constant in space and time whereas here the re- 
sistivities locally increase significantly during the 
course of the simulations. A reasonable inference 
is that the influence of local variations in resis- 
tivities averages out in some sense on the global 
scale. 

As we shall see later, however, the effect of spa- 
tially varying resistivities is noticeable in proper- 
ties such as the power spectrum of the density and 
magnetic field. 

The various energy decay rates can be modeled 
approximately as power-laws in time, i.e. t~@ . Fit- 
ting the kinetic energy, /3k , magnetic energy, /3b, 
and total energy, /3xot , as functions of time in this 
way we obtain the values given in table [2] This 
data confirms quantitatively what can be observed 
in Figure [5] and extends it to the decay of the en- 
ergy in magnetic perturbations: increasing reso- 
lution decreases the rate of energy decay, but the 
difference between the 256 3 and 512 3 simulations 
is relatively minor. We note in passing that the de- 
cay in the energy in magnetic perturbations is con- 
siderably more sensitive to resolution than kinetic 
energy: /3b varies between 1.43 and 1.28 while /3k 
only varies between 1.38 and 1.34. We know that 
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Fig. 2. — Plot of the normalized ambipolar re- 
sistivity (top panel), normalized Hall resistivity Fig. 3. — As figure [5] but the snapshots are taken 
(middle panel) and r) = t^t. Note that the scale at t — 0.2 i c . Note that the scale here is logarith- 
is logarithmic in the top and middle panels while mic in each case including the bottom panel, 
it is linear for the bottom panel. 
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the resistivities arc a critical factor in determining 
the decay of the magnetic energy since they facili- 
tate loss of magnetic energy through reconnection. 
We attribute the extra sensitivity to resolution of 
/3b to the necessity to properly resolve the diffu- 
sive effects in the induction equation, including the 
variation of the resistivities themselves. The vari- 
ation of the resistivities throughout the computa- 
tional domain is significant (see figure [2]) and it 
is interesting to note that /3b in simulations with 
fixed resistivities is not so sensitive to resolution 
(see Paper I). 

4.2. Energy decay 

We now discuss the behavior of the kinetic and 
magnetic energy in our multifluid simulations and 
compare with those in Paper I. 

4- 2.1. Kinetic energy decay 

Figure [5] contains plots of the decay of kinetic 
energy with time for the 512 3 resolution simula- 
tions outlined in table [T] Note that the kinetic en- 
ergy decay in all of the simulations is very similar 
until around t rj 0.02 t c . This is because at such 
early times compressions are only just starting to 
form and so the non-ideal terms in the induction 
equation have had almost no effect on the dynam- 
ics. The subsequent energy decay of simulations 
mhd-512 and hall-512 arc almost identical to each 
other. The energy decay of simulations mc-512 
and ambi-512 are virtually identical over the full 
plotted range while the data plotted for fr-512 co- 
incide with the former simulations for times in the 
range 0.03 t c < t < 0.2t c . We have fitted the ki- 
netic energy decay as a power-law in the range 
0.2t c < t < t c (i.e. after approximately one initial 
flow crossing time) and the exponents are given in 
the first column of table [2] 

It is clear that the presence of ambipolar dif- 
fusion has a significant impact on the behavior of 
the kinetic energy in the turbulent system. This 
is a result of the exchange of energy between ki- 
netic and magnetic energies as will be discussed in 
section K2~2l 

From figure [5] and table [2] it is evident that the 
Hall effect has almost no impact on the kinetic 
energy decay in turbulence in molecular clouds. 
In order to emphasize any possible impact of the 
Hall effect we have plotted the time evolution of 
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Fig. 4. — Log scale plot of the kinetic energy (nor- 
malized to its initial value) for each of the simula- 
tions in the resolution study. 

Table 2: The values of the exponent for the kinetic, 
magnetic and total energy decay for the simula- 
tions presented in this work. These exponents are 
calculated by fitting the data over the time inter- 
val [0.2t c ,* c ]. 



Simulation 


Pk 


/3b 




mc-64 


1.38 


1.43 


1.39 


mc-128 


1.37 


1.35 


1.36 


mc-256 


1.35 


1.30 


1.33 


mc-512 


1.34 


1.28 


1.32 


hall-512 


1.12 


1.05 


1.09 


ambi-512 


1.35 


1.30 


1.33 


mhd-512 


1.12 


1.06 


1.10 


fr-512 


1.30 


1.28 


1.30 




0-001 1 1 1 1 — 1 1 1 — 1 

0.01 0.1 1 

Time (crossing limes) 

Fig. 5. — Log scale plot of the kinetic energy (nor- 
malized to its initial value) for each of the 512 3 
simulations. 
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the ratio of the kinetic energy in each of our sim- 
ulations to that in mc-5-512 in figure[6]on a linear 
scale. It is clear even in this figure that the Hall 
effect has little impact on the evolution of the tur- 
bulence. This result is also reproduced if we ex- 
amine the evolution of the magnetic energy (see 
section 14.2.2)) . This supports our conclusion from 
Paper I in which the simulations were run using 
fixed resistivities (see also fr-512 in this work). 

Also shown in figure [5] is the energy decay for 
simulation fr-512. Given the wide variation of the 
resistivities in both space and time (see figures 
[2] and [3} it is somewhat surprising that the en- 
ergy decay is so similar to that of mc-512. In fact, 
the volume average of the ambipolar resistivity at 
t = 0.5 t c in mc-512 is approximately 60% higher 
than that in the fr-512 simulation. It would appear 
that, while ambipolar diffusion enhances energy 
loss, the expected spatial and temporal variation 
of it does not have that much influence. 

One possible reason for the small divergence be- 
tween mc-512 and fr-512 would be if the locations 
in which the resistivity is high are regions in which 
the magnetic field, B, is varying weakly. To ex- 
plore this we define a scalar, SB, by 

SB ee |grad(£y| + |grad(B„)| + |grad(B 2 )| (25) 

where grad^^), for example, is a normalized gra- 
dient defined by 

, / S X B X \ 
grad^) ee — S y B x (26) 
Bo V S Z B X J 

where S. means centered differencing in the indi- 
cated direction without normalizing by the zone 
spacing and Bq is the magnitude of the magnetic 
field throughout the domain at t = 0. SB is then 
a dimensionless measure of the variation of B at 
any point in space or time: if SB is large it means 
there is a high gradient in one or more of the com- 
ponents of B and therefore resistivity will have an 
important influence here. Figure [7] contains snap- 
shots of SB and |B| at t = 0.5 t c . There is rich 
structure in SB which is not apparent in |B|. 

Figure[5]contains a histogram plot of the two di- 
mensional probability density function for SB and 
ta (the ambipolar resistivity). What is striking 
about this plot is that there is a notable lack of 
high SB with corresponding high resistivity. Of 
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Fig. 6. — Plot of the ratio of the kinetic energy in 
our 512 3 simulations to that in mc-5-512 




Fig. 7. — Slices showing the z = 0.5 plane at 
t = 0.5 t c for both |B| (left panel) and SB (right 
panel) . 
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course, if we have any system in which there are 
regions of high and low resistivity we expect that, 
over time, the regions with high diffusion will have 
lower variation in B so this, in itself, doesn't tell 
us much. However, it does prompt us to look a 
little more closely at the behavior of r\. 

Comparison of the middle panel of figure Q] with 
the top panel of figure [5] shows that the ambipo- 
lar resistivity is higher in regions of low density, 
as would be expected. Now, in supersonic/super- 
Alfvenic turbulence kinetic energy is dissipated 
most strongly at strong shocks. However, shocks 
propagating into regions of very low density will 
not dissipate kinetic energy effectively. Since it is 
precisely these regions in which our resistivities are 
high we must conclude that, in fact, the regions 
of enhanced resistivity do not contribute signifi- 
cantly to energy dissipation and hence we would 
not expect the introduction of spatially varying 
resistivities to increase the rate of energy decay 
in our simulations. Further, since the volume av- 
erage of, for example, the ambipolar resistivity is 
actually higher than that used in fr-512 we would 
not expect the spatial variation of this resistivity 
to reduce the rate of energy decay cither. 

4-2.2. Magnetic energy decay 

We now move on to discuss the decay in mag- 
netic energy. Initially, as outlined in Paper I, the 
magnetic energy increases as the flow compresses 
and stretches the magnetic field throughout the 
computational domain. Once this initial increase 
in the energy has occurred it is gradually lost 
through two main avenues: magnetic reconnection 
and transfer of magnetic energy to kinetic energy 
which can then be dissipated in shocks and other 
viscous processes. 

Figure |9] contains plots of the decay of mag- 
netic energy with time. It is clear that, in common 
with the case of kinetic energy, the hall-512 and 
mhd-512 simulations arc almost identical while 
the ambi-512 and mc-512 simulations arc also well 
matched. This supports our inference from sec- 
tion 14.2.11 that the Hall effect has almost no im- 
pact on energy decay in molecular clouds on the 
global scale. 

Ambipolar diffusion does, however, have a sig- 
nificant impact on the behavior of both the kinetic 
and magnetic energies. This was also noted in Pa- 




Fig. 8. — Probability density function for the nor- 
malized variation in B and the resistivity normal- 
ized to its initial value (i.e. the same as that used 
in simulation fr-512). It is clear that high resistiv- 
ities only occur where there is low variation in B. 
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Fig. 9. — As figure [5] but for magnetic energy. 
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per I and we explain it in the same way, recapitu- 
lated briefly here for completeness. Consider a re- 
gion of the flow undergoing compression. During 
this compression kinetic energy will be converted 
into both magnetic and internal energy through 
increasing the magnetic pressure and the thermal 
pressure. Given that this is a turbulent flow we ex- 
pect that, after some time, this region will begin 
to expand again. However, during the compres- 
sion ambipolar diffusion will have diffused away 
some of the magnetic energy thereby leaving less 
to be converted back to kinetic energy. In this way 
the presence of ambipolar diffusion creates a new 
path through which energy can be lost from the 
system and reduces the level of all forms of energy 
in the system. 

4.3. Power spectra 

We now move on to a study of the power spec- 
tra obtained from the multifluid MHD turbulence 
simulations. These spectra are important from the 
point of view of understanding the types of struc- 
tures formed by the turbulence and are a more 
discerning tool for exploring any structural differ- 
ences caused by multifluid effects. Table [3] con- 
tains the exponents of the power spectra assuming 
a power-law relationship between power and wave 
number. All the analysis presented in this section 
is performed on data taken at t = t c . 

4-3.1. Density power spectra 

We turn first to the scales of the structures 
formed in the density distributions for our vari- 
ous simulations. Figure [10] contains plots of the 
power spectra of the neutral density (or density, 
in the case of simulations mhd-512 and fr-512) for 
all the 512 3 resolution simulations. 

The power spectra for the simulations includ- 
ing the effects of ambipolar diffusion are approx- 
imately broken power laws made up of 3 distinct 
power laws: l<fc<10, 10<fc<100 and 100 <k. 
Below the low k break, the spectrum is dependent 
on the scale of the computational domain. At high 
k approaching the grid scale, numerical viscosity 
will begin to dominate. In common with the re- 
sults presented so far we see that there is little 
difference between simulations mc-512 and ambi- 
512 - in fact it is difficult to distinguish between 
the two spectra without careful examination of fig- 



ure [TO] The fr-512 power spectrum is also similar 
to mc-512 and ambi-512 although it has slightly 
less power at intermediate values of k with the 
difference here being at most 10%. 

There is almost no detectable difference be- 
tween simulations hall-512 and mhd-512. Evi- 
dently, the Hall effect has a much weaker influence 
on the density structure in molecular clouds than 
ambipolar diffusion. Ambipolar diffusion, on the 
other hand, has a very significant impact with pro- 
nounced damping of density structures at scales 
less than one tenth of the domain size (correspond- 
ing to a physical scale of approximately 0.02 pc). 
This damping is evident in figure [JJ where the den- 
sity structures in simulations with ambipolar dif- 
fusion are more smeared than in mhd-512. 

From a quantitative perspective, the data in ta- 
ble [3] shows that the inclusion of ambipolar diffu- 
sion significantly steepens the power spectra, in- 
creasing the exponent by more than 0.5 over the 
cases which do not include the effect. The data 
again suggests that the Hall effect has minimal 
impact. 

The results for fr-512 show a softer spectrum 
at large length scales and a harder spectrum at 
short scales than mc-512 and ambi-512, indicat- 
ing that self-consistent calculation of the resistivi- 
ties reduces the level of fine structure. This result 
must, however, be confirmed by higher resolution 
simulations before it can be regarded as reliable. 

Figure [IT] contains plots of the power spectra 
of the neutral density and the density of the neg- 
atively charged species at t = t c for comparison. 
It can be seen that the neutral mass density has 
more power for k > 10 although the qualitative 
shape of the power spectra are the same in each 
case. This is in q ualitative agree ment with the re- 
sults presented in lLi et al.l (|2008l ) for driven turbu- 
lence simulations. For comparison with the results 
in Table [3] the exponent for the charged species 
mass density in the range 4 < k < 10 is -2.21 
while in the range 10 < k < 100 it is -4.84. The 
power spectrum for the neutral mass density is 
harder than that for the magnetic field and, in the 
range 10 < k < 100, the spectrum for the charged 
species mass density lies somewhere between the 
two. This is not surprising as it is a function of the 
neutral density and velocity (through drag) and 
the magnetic field (through the Lorentz force). 
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Table 3: The values of the exponent for the power spectra of density, velocity and magnetic field measured 
at t = t c . All fits are over the range 5 < k < 20 unless otherwise noted. 



Simulation 


Density 


Velocity 


Magnetic field 


mc-512 


1.82 a , 4.33 b 


1.34 


2.14 a , 5.38 b 


ambi-512 


1.79 a , 4.31 b 


1.34 


2.15 a , 5.43 b 


hall-512 


1.25 


1.01 


1.51 


mhd-512 


1.27 


1.00 


1.55 


fr-512 


2.11 a , 4.21 b 


1.28 


2.20 a , 5.53 b 



"Fitted over 4 < k < 10 
^Fitted over 10 < k < 100 




Fig. 10. — Plot of the normalized, spherically inte- 
grated power spectrum of the neutral mass density 
for each of the 512 3 simulations at t = t c . 



Charged Density - 
Neutral Density - 




Fig. 11. — Plot of the normalized, spherically in- 
tegrated power spectrum of the neutral and neg- 
atively charged mass densities for simulation mc- 
512. 



4-3.2. Velocity power spectra 

Figure [T2] contains plots of the velocity power 
spectrum for the neutral velocity. In common with 
section 14.3.11 those simulations incorporating am- 
bipolar diffusion have significantly less power at 
almost all scales than those without. This is to be 
expected given the increased rate of loss of turbu- 
lent energy in the presence of ambipolar diffusion 
(see section |4~2")) . 

As found in Paper I, and again here, there are 
clear differences of a qualitative nature between 
the density power spectra and the velocity power 
spectra for the multifluid simulations. Those sim- 
ulations incorporating ambipolar diffusion (mc- 
512 and ambi-512) exhibit a strong power-law in 
the range 10 < k < 100 in the density power 
spectrum. This is not true of the velocity power 
spectra. For the latter spectra there is a break at 
roughly k ~ 20 and again at k ~ 100. The lat- 
ter break we interpret as being at the scale where 
numerical diffusive effects begin to dominate the 
non-ideal effects in the induction equation. The 
lower break can reasonably be interpreted as the 
scale at which the non-ideal effects become impor- 
tant. The marked qualitative differences between 
the density and velocity power spectra indicate 
a considerable decoupling between the two vari- 
ables. It is worth recalling here that the power 
spectra are calculated at t = t c so it is reasonable 
to expect that the turbulence is well developed at 
this stage. 

In common with the density power spectra, the 
presence of ambipolar diffusion produces steeper 
velocity power spectra (see T able EH in qualit ative 
agreement with the results of iLi et all ( 2008h . 
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The results at high k are, of course, dominated 
by numerical diffusive effects. However, it is in- 
teresting to note that simulation fr-512 actually 
has very slightly less power at short scales than 
mc-512 despite the volume average of ta being 
roughly twice as high as the resistivity used in fr- 
512. This adds weight to the inference from sec- 
tion 14.2.11 that regions of high resistivity tend not 
to be coincident with regions of high gradients in 
the magnetic field and therefore do not have the 
level of influence one would naively expect. 

FigurclT3"l contains plots of the power spectra for 
the velocity of the neutral and negatively charged 
species. It can be seen that, except at very high 
k 7 they are virtually identical. Given that the 
charged velocity is defined by balance between 
drag with the neutrals and the Lorentz force this 
is, perhaps, unsurprising. 

4- 3. 3. Magnetic field power spectra 

Figure fT?] contains plots of the spherically in- 
tegrated power spectra for the magnetic field. 
Once again, ambipolar diffusion has a much big- 
ger impact on these spectra than the Hall effect. 
The magnetic field power spectra are considerably 
steeper at all k in its presence. The absolute power 
at any scale is also considerably lowered by am- 
bipolar diffusion, as would be expected from the 
discussion in section 14.2.21 There is a qualitative 
similarity between the density power spectra (fig- 
ure llO[) and the magnetic field power spectra which 
is absent when comparing the latter with the ve- 
locity power spectra (see figure [1~2")) . 

Interestingly, we see the phenomenon that fr- 
512 has very slightly less power at k > 20 than 
mc-512 indicating that fixing the resistivities at 
t = actually results in slightly more dissipation 
than allowing it to vary in time and space. 

We find that the Hall effect has a slightly more 
noticeable effect on the magnetic field power spec- 
tra than on the spectra of velocity or density: it 
gives rise to a little more structure on short scales 
which is absent in its absence. This would be ex- 
pected as the Hall effect is a dispersive effect act- 
ing directly on the magnetic field. The results for 
the density and velocity power spectra, however, 
demonstrate that this influence over the magnetic 
power spectra does not translate into an influence 
over the other variables. 




Fig. 12. — Plot of the normalized, spherically inte- 
grated power spectrum of the neutral velocity for 
each of the 512 3 simulations at t = t c . 



CharsjL'd Velocity 
Neutral Velocity 



Fig. 13. — Plot of the normalized, spherically inte- 
grated power spectrum of the neutral velocity and 
the velocity of the negatively charged species for 
simulation mc-512 at t = t c . 
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4.4. Velocity dispersion 



It has been widely reported that the observed 
velocity dispersion in molecular clouds behaves as 
a power-law in the size of the field of view. Figure 
[15]contains plots of the velocity dispersion for each 
of the 512 3 resolution simulatio ns presented in this 
work As noted in Paper I, and lLemaster fc Stone 
(l2009h . no power-law is observed. This may be 
due to the fact that in (multifluid) MHD turbu- 
lence there are many relevant signal speeds due 
to the variation in the orientation and intensity 
of the magnetic field, in contrast to the situation 
with hydrodynamic turbulence. Hence there is no 
reason to expect to see a power-law. Finally, it 
should be noted that, at this point in the simu- 
lation, the turbulence has decayed such that it is 
largely sub-sonic or transonic. This may also ex- 
plain the lack of a power-law. Indeed, it should 
also be noted that our results here are not neces- 
sarily in contradiction with observations since the 
velocity dispersion-size correlation can only be ac- 
curately measured in the supersonic regime. 

It is clear that ambipolar diffusion reduces the 
velocity dispersion at all length scales - this is 
what would be expected given the results from the 
velocity power spectra presented in section 14.3.21 
Once again the results here indicate that the Hall 
effect has little impact. Again, it is interesting to 
note that the spatial variation of the resistivity in 
mc-512 appears to have almost no impact as the 
velocity dispersion seen in fr-512 is almost identi- 
cal to that in mc-512. 

5. Conclusions 

We have presented results from a suite of 512 3 
resolution simulations of fully multifluid MHD de- 
caying turbulence. The effects incorporated in- 
clude the Hall effect and ambipolar diffusion. We 
have performed a resolution study to ensure that 
the energy decay rate, being the main result pre- 
sented here, is reliable. We have confirmed the re- 
sults of the simplified calculations in Paper I that 
the Hall effect has little impact on the nature and 
behavior of turbulence in molecular clouds under 
the well motivated physical parameters assumed 
in this work. Further, the presence of ambipo- 
lar diffusion increases the rate of energy decay at 
length scales of 0.2 pc and less. The same conclu- 
sion is drawn for the behavior of the energy in the 




Fig. 14. — As figure [12] but for the magnetic field. 




Fig. 15. — Log-scale plot of the velocity dispersion 
as a function of length scale for each of the 512 3 
simulations. There is no power-law apparent in 
this data. 



14 



magnetic field. 

The power spectra for these simulations again 
suggest that the Hall effect has little impact on the 
flows with the exception of the spectrum of mag- 
netic field variations. We must keep in mind that 
the maximum resolution used here (512 3 ) is only 
enough to resolve about half a decade in fc-space 
and it is therefore difficult to be confident of the 
details of the power spectra. Notwithstanding this 
consideration it does appear clear that ambipolar 
diffusion steepens the power spectra of the neu- 
tral velocity, density and the magnetic field. As 
noted in Paper I, it appears that at a resolution 
of 512 3 and an assumed length scale of 0.2 pc we 
have resolved the length at which ambipolar dif- 
fusion begins to influence the turbulent cascade. 
In Paper I only constant resistivities were imple- 
mented and hence it was unclear whether this lat- 
ter result would survive the inclusion of more re- 
alistic fully multifluid MHD in which the resistiv- 
ities vary strongly in space and time. The results 
presented here imply that it does. 

The power spectra of the neutral velocity and 
the magnetic field differ qualitatively from that of 
the density with breaks occurring in the former 
which are not seen in the latter. This suggests a 
decoupling between these fields. 

The velocity dispersion as a function of length 
docs not behave as a power law. This is not unex- 
pected as the nature of MHD turbulence implies a 
wide range of applicable signal speeds which can, 
when combined, remove the power-law behavior 
which might be expected if only one signal speed 
were relevant. 
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